Extraordinary transmission caused by symmetry breaking 
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The terahertz transmission properties of symmetrical and asymmetrical annular apertures arrays 
(AAAs) are investigated both experimentally and numerically. It is found that only odd order 
resonant modes are observed for the symmetrical structures but both odd and even order resonances 
can be shown for the asymmetrical structures. Breaking the symmetry of AAAs by gradually 
displacing the H-shaped AAAs to U-shaped AAAs allows an intensity modulation depth of 99% 
of the second order resonance. Simulation results verify the experimental conclusions well. This 
result provides a tremendous opportunity for terahertz wavelength tunable filtering, sensing, and 
near-field imaging. 
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Since the extraordinary transmission (EOT) phe- 
nomenon of sub- wavelength metallic apertures arrays has 
been reported by Ebbesen and his co-workers [H, 0, it 
has been paid extensive attentions to understand the un- 
derlying physical mechanisms 0, Q ■ Although it is still 
debated on the real mechanism of the EOT, surface plas- 
mon polaritons (SPPs) originating from the coupling of 
light with the surface charges oscillation on the metal- 
dielectric interface Q and the localized surface plasmon 
resonance (LSPR) excited on the hole ridge are recog- 
nized as the main possible contributions to the EOT [J] . 
Moreover, it has been shown that the transmission spec- 
tral response of sub- wavelength metallic apertures arrays 
are affected by the cell size Q, cell shape d, Q, lat- 
tice periodicity @, metal film thickness @, the refrac- 
tive index of the adjacent medium as well as the 
polarization of incident light [Tl| . Annular aperture ar- 
rays (AAAs) have also been demonstrated to have EOT 
(la Il3j . The influence of geometric and structure pa- 
rameters of the AAAs on the EOT properties has been 
extensively investigated. However, the influences of the 
symmetry of the AAAs on the EOT have not been re- 
ported. In this Letter, the transmission properties of 
the symmetric and asymmetric AAAs are investigated 
both experimentally and theoretically. The interactions 
of the terahertz (THz) electromagnetic radiation with the 
AAAs have been measured with the traditional THz time 
domain spectroscopy (TDS). The field distributions of 
the resonant wavelengths are simulated with the finite 
difference time domain (FDTD) method. The physical 
original of the resonant enhanced transmission has been 
analyzed. It is found that the symmetry breaking can 
enhance the even order mode EOT. 
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FIG. 1: Schematic view of the H-shaped AAAs. (a) Top 
view and (b) side view of the sample. The x-polarized THz 
illumination is also indicated in the figure. 



The proposed AAAs structure is schematically shown 
in Fig. [TJ It is a H-shaped annular slit. The outer arm 
length of the H-shaped slit is L\, the upper and lower in- 
ner arm lengths are L2 and L3, the outer and inter arm 
widthes are W\ and W2, respectively. The slit width is 
W3 — 2.5fj,m and the period (A) of the sample is 200/^m. 
This sample is symmetrical along both x and y directions. 
The total average circumference of the AAAs, which is 
the average length of out and inner edge of the slit, is 
denoted as L to tai- The samples are fabricated onto a 
thin copper film (t — 200nm) with the conventional pho- 
tolithography and metallization processes. The substrate 
is a 500/xm thick doubly polished p-type silicon wafer. 

The transmission spectra of the samples are measured 
by using a typical THz TDS. The THz wave pulses are 
generated from a dipole-type photoconductive antenna 
illuminated by a 100/s laser with 800nm center wave- 
length. The valid range of the system is 0.2 — 2.QTHz. 
The detection is achieved by a ZnTe crystal. For the ref- 
erence, a bare silicon wafer identical to the one on which 
the AAAs are fabricated is used. The transmitted THz 
pulses are measured at normal incidence such that the 
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FIG. 2: Measured (symbol) and simulated (solid line) trans- 
mission spectra for AAAs with same average circumference 
Ltotai but different L 2 . 



electric and magnetic field of the incident radiation arc 
in AAAs' plane. All measurements are done at room 
temperature and in a dry atmosphere to mitigate water 
absorption. The normalized transmittance is obtained 
by T(y) = \E samp i e {v) / E ref (v)\, where E sample (y) and 
E re f (v) are Fourier transformed transmitted electric field 
of the sample and reference pulses, respectively. 

The transmission spectra of five samples with different 
Li are shown in Fig. [5] with symbols. The parameters 
of the structures are selected as follows: The widthes W\ 
and W% are 30/im and 20//m, respectively. The length L\ 
is 80/im. The length Li increases from 20fj,m to 40/im 
with step of 5[im and thus Lz reduces from 20/im to 
Ofim. The average circumference of the slit is kept as a 
constant as L to tai = 390/im. For Li = L% = 20/xm, the 
sample is H-shaped AAAs and symmetrical along both 
x and y directions. For L 2 = 40/zm and L 3 = Ofim, the 
sample becomes to U-shaped AAAs and it is only sym- 
metry along the x direction. It can be seen that two 
resonance peaks at 0.35THz and l.HTHz appear in all 
transmission spectra. With Li increasing, the symmetry 
of the structure along the y direction has been broken, 
the locations of these two peaks nearly unchanged but a 
new resonance peak appears at 0.70THz. The transmis- 
sion for 0.7QTHz is enhanced from 1.97% for symmetri- 
cal structure to 17.5% for asymmetrical one, the intensity 
modulation depth is 99%. 

For the AAAs structure, the localized surface plas- 
monics will be excited under the light illumination. The 
surface charges are assembled at the edge of the metal. 
Based on the basic idea of Ref. [lij], the resonance ap- 
pears when the following condition is satisfied: 
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where f m is the resonant frequency, m is an integer which 
is the order of the resonant mode, c is the speed of light in 
vacuum, L to tai is the average circumference of the AAAs, 
and n e f / is the effective refraction index for the localized 



surface plasmonic wave. In our case, n e ff is approxima- 
tively equal to (1.0 + n. su &)/2, n su b is the refraction index 
of the substrate. Under this condition, the standing wave 
of the localized surface plasmon can be generated along 
the wedges of the metallic film. Furthermore, the dis- 
tribution of localized surface plasmon should also satisfy 
the symmetry requirement of the samples. All samples 
considered here are symmetrical along the x direction 
and the illumination light is polarized along the x di- 
rection, thus the distribution of the assembled surface 
charges should be antisymmetrical along the x direction. 
If the sample is also symmetrical along the y direction, 
the distribution of the assembled surface charges should 
also be even symmetrical along the y direction. When 
the distribution of the assembled surface charges satisfy 
both the standing wave condition and symmetry condi- 
tion, the corresponding transmission can be greatly en- 
hanced. For all samples considered here, the odd mode 
(for example, m = 1 or m = 3) can easily satisfies these 
requirements no matter the sample is symmetrical along 
the y direction or not. However, if the resonant mode 
is even, the symmetry of the sample along the y direc- 
tion will affect the distribution of the surface charges and 
hence the transmission amplitude. It requires the distri- 
bution of the surface charges is also symmetrical along 
the y direction. Thus, the standing wave condition for 
second mode can not be stratified any more, therefore, 
the transmission for the second resonant mode will be 
greatly reduced, as shown in Fig. 2 with Li — 20fim. 
But when the symmetry along the y direction is broken, 
the requirement of the symmetry of surface charge distri- 
bution along this direction does not exist, the standing 
wave of the localized surface plasmon can be generated 
and the resonant peak appears, as shown in Fig. [2] with 
Li = 40/im. 

In order to get an insight into the nature of the res- 
onant transmission for the proposed structures, the nu- 
merical simulations are carried out by using the FDTD 
approach [l5[ . The metal copper is described by the per- 
fect electric conductor to fit its realistic characteristic in 
the terahertz range. The refraction indices of the air and 
substrate silicon are set as 1.0 and 3.31 in the simula- 
tions, respectively. The calculated transmission spectra 
are also shown in Fig. [5] with solid lines. It can be found 
that the simulation results correspond to the experiment 
measurements well. From the simulation results, it can 
also be found that the sample with Li = 20/im has two 
resonant peaks at Q.35THz and l.HTHz, which corre- 
spond to m = 1 and m = 3, respectively. But the sample 
with Li — 40/im has three resonant peaks. A new peak 
corresponding to m = 2 appears at 0.70THz. With Li 
increasing, the transmission for the second resonant mode 
enhanced. 

Figure [3] presents the extremum electric field distribu- 
tions at the aforementioned three resonant frequencies 
for samples with Li equals 20/um, 30/zm, and 40/im, re- 
spectively, which exhibits the distributions of the surface 
charges. When the structure is symmetrical along both 
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FIG. 3: Simulated E distributions on the interface between 
metal film and dielectric substrate silicon for three samples at 
three resonant wavelengths. The snapshot moment is when 
the field amplitude reaches its maximum. Only a unit cell on 
the x-y plane is shown. The field quantity is normalized with 
respect to the incident field Eq. 



x and y directions, it can be found from Fig. |3ta)-(c) 
that the extremum electric field distributions are anti- 
symmetrical along the x direction but symmetrical along 
the y direction. These are natural results for satisfying 
the symmetrical requirement of the samples. For the odd 
modes m — 1 and m = 3 (correspond to Fig. 02a) and 
(c), respectively), the standing waves are generated along 
the edges of the metallic film. One wave crest along the 
out side of slit in Fig. [UJa) and three crests in Fig. [3jc) 
can be found. However, in order to satisfy the symmetry 
condition along the y direction, the surface charge wave 
can not generate a standing wave with mode to = 2, as 
shown in Fig. 3(b). Only a temporarily electric field dis- 
tribution is shown here and no standing wave distribution 
can be found. For the structure with L 2 = 40/xto, the 
electric field distributions are also antisymmetrical along 
the x direction. Since the structure is not symmetrical 
along the y direction any more, the surface charges dis- 
tributions are not required to be symmetrical along the 
y direction. There are 1,2, and 3 standing waves gener- 
ated for mode to = 1 (Fig. 3(g)), m = 2 (Fig. 3(h)) and 
to = 3 (Fig. 3(i)), respectively. The even mode standing 
wave is formed and thus the corresponding transmission 
has been enhanced. It can be concluded that the reso- 
nant peak to = 2 for the structure with L2 = 40^tm is 



caused by the symmetry breaking along the y direction. 

For further validation, the transmission spectra for an- 
other series of structures are calculated and shown in Fig. 
[U The AAAs is cross shaped and the average circumfer- 
ence L c totai is still 390/im. The girder is shifted from the 
middle to the bottom of the cross to break the symmetry 
along the y direction. It can be seen from the transmis- 
sion spectra that the odd mode m = 1 appears around 
0.35THz for all structures. However, the transmission 
intensity for second mode to = 2 increases from 0.39% 
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FIG. 4: (a) Schematic of cross shaped AAAs with fixed pa- 
rameters: Wd = Wd = 100/im, Wcz = 2.5/im, L c i + 
L C 2 = 80/im, L C 3 = 20/im, Lctotai = 390/xm, and the period 
A = 200/tm. (b) Simulated spectra for cross shaped AAAs 
with with same average circumference L cto t a i but different 

L C 2- 



for symmetrical structure to 6.5% for the asymmetrical 
structure. These results also demonstrate that the sym- 
metry breaking can enhance the odd mode resonance. 

In summary, the transmission characteristics of the 
symmetrical and asymmetrical AAAs with narrow gap 
are investigated experimentally and numerically. The 
measured and simulated transmission spectra reveal that 
the resonant transmission are affected by the geometric 
structure. Both even and odd order modes can emerge 
for the asymmetrical AAAs but only odd order modes 
appears for the symmetrical AAAs. The appearance of 
the even mode is caused by the symmetry breaking. This 
research indicates a new freedom in tuning the electro- 
magnetic response, which offers a path to design more 
robust plasmonic devices. 
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